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Synchronization of networked oscillators is known to depend fundamentally on the interplay be¬ 
tween the dynamics of the graph’s units and the microscopic arrangement of the network’s structure. 

For non identical elements, the lack of quantitative tools has hampered so far a systematic study of 
the mechanisms behind such a collective behavior. We here propose an effective network whose topo¬ 
logical properties reflect the interplay between the topology and dynamics of the original network. 

On that basis, we are able to introduce the ’’synchronization centrality”, a measure which quantifies 
the role and importance of each network’s node in the synchronization process. In particular, we 
use such a measure to assess the propensity of a graph to synchronize explosively, thus indicating 
a unified framework for most of the different models proposed so far for such an irreversible transi¬ 
tion. Taking advantage of the predicting power of this measure, we furthermore discuss a strategy 
to induce the explosive behavior in a generic network, by acting only upon a small fraction of its 
nodes. 

PACS: 89.75.Hc, 05.45.Xt, 87.18.Sn, 89.75.-k 

One of the most intriguing processes in complex net¬ 
works’ dynamics is synchronization: the spontaneous or¬ 
ganization of the network’s units into a collective dy¬ 
namics. This phenomenon is known to be related to a 
delicate interplay between the topological attributes of 
the network and the main features of the dynamics of 
each graph’s unit m • The conditions for synchroniza¬ 
tion in complex networks have been addressed by means 
of different approaches. For identical units, one of the 
most successful tools is, for instance, the Master Stabil¬ 
ity Function [4], which rigorously shows how the spectral 
properties of the graph influence the stability of synchro¬ 
nization [I]. However, the general case of non-identical 
units is far more complicated, and often needs a numer¬ 
ical approach, where the topology-dynamics relationship 
can only be investigated within specific scenarios W- 

Such a connection between structure and dynamics of 
a network is of particular importance in the case of the 
recently reported explosive synchronization (ES), an irre¬ 
versible and discontinuous transition to the graph’s syn¬ 
chronous state. Originally, ES was described in all-to- 
all coupled ensembles of Kuramoto oscillators jl.0] for a 
specific distribution of natural frequencies D-p. Later 
on, various kinds of degree-frequency correlations were 
found to be able to induce ES in networks of periodic 
and chaotic oscillators EUHE], or neural networks [IB] . In 
addition, other microscopic mechanisms were proposed, 
based on diverse coupling strategies P2H31], or by intro¬ 
ducing adaptive dynamics in a fraction of the network’s 
units [20]. 

In this Letter, we propose the use of an effective net¬ 
work whose structure indeed reflects the interplay be¬ 
tween the topology and dynamics of the original system. 

On that basis, we introduce the synchronization central¬ 
ity as a measure to quantify the role of each node in the 


synchronization process. This measure allows us to pro¬ 
vide a general explanation of the mechanisms underlying 
ES and revisit the main scenarios where such a behavior 
was previously reported. Finally, we formulate a crite¬ 
rion to induce explosive transitions by acting only on a 
small fraction of the network’s nodes. 

We start by considering a network of N phase oscilla¬ 
tors, whose instantaneous phases evolve in time according 
to the Kuramoto model m : 

N 

Oi = Ui + ip) E Aij sin (0j - 0i) i = l,-,N, (1) 

where Oi is the phase of the i-th oscillator, Ui its natu¬ 
ral frequency (chosen from a generic, known, distribution 
a the coupling strength, and (.) stands for the en¬ 
semble average. The topology of the network is encoded 
in the adjacency matrix A (A^ = 1 if node i is linked to 
node j, and = 0 otherwise). The degree of a node is 
ki = Yj Aij . The level of synchronization is measured 

by the order parameter r = '4f(\YiLi eGi \)T, where |.| 
and (.)t denote module and time average, respectively. 
Along this Letter, the network size is fixed to N = 1, 000 
and the natural frequencies uoi are randomly drawn from 
a uniform distribution in the interval [—0.5, 0.5], unless 
otherwise specified. 

As the coupling strength a gradually increases, system 
0 experiences a transition from an incoherent (r ~ 0) 
to a synchronized state (r ~ 1), a process often referred 
to as path to synchrony (PTS) [21]. In heterogeneous 
networks, the PTS is mainly dominated by the most con¬ 
nected nodes (or hubs) which actually act as synchroniza¬ 
tion seeds, and progressively recruit the other network’s 
nodes. At variance, in homogeneous networks, the PTS 
is characterized by the emergence of coherent clusters 
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growing around multiple synchronization seeds. Along 
the PTS, the first nodes that locally synchronize gener¬ 
ally correspond to pairs of connected oscillators whose 
natural frequencies are closer, whereas the globally syn¬ 
chronized state emerges around those with natural fre¬ 
quencies close to the synchronization frequency fl s . 

While traditionally attention has been focused on the 
natural dynamics of each node, we now discuss a dif¬ 
ferent approach, where the links prevail over the nodes 
themselves, and we use the frequency detuning A ujij = 
\uji — ujj | between each pair of nodes as the key dynam¬ 
ical feature for the determination of the PTS. To for¬ 
malize our idea, let us introduce a change of variables 
r i e l ^ i = Xqeiy ? where Vi(t) is a local order pa¬ 
rameter and Ti is the set of neighbors of the node i. 
Substituting into Eq. 0 we obtain 

6i = uji + ari sin (\Eq — Of) . (2) 



FIG. 1. (Color online). Synchronization transition in per¬ 
turbed (a) ER, < k >= 30 and (b) SF < k >= 12 networks, 
as compared with the original unperturbed system (black cir¬ 
cles). Red triangles and blue squares shows that a positive 
(negative) perturbation of the adjacency matrix with S = 0.5 
(S = —0.5) reinforces (weakens) the local coupling of those 
links with large detuning, and results in an increase (frustra¬ 
tion) of the level of synchronization. Furthermore, the local 
perturbation enhances/frustrates the synchronization more 
efficiently than a global perturbation over all links (dashed 
lines, where SAuJij —> 5 {Aw), being (Au) the average over 
nonzero values of the detuning matrix). 


It follows naturally that the velocity difference is = 
0i — 0j = Wi—Wj+a [ri sin(\Jq — Of) — rj sin(Tj — 0j)\. As 
synchronization implies = 0, the set of links through 
which synchronization may take place must fulfill 


Auij < cr (n+Tj), (3) 

which in fact relates local synchronization to the fre¬ 
quency detuning associated to the links, rather than 
to the microscopic structure of connections around the 
nodes. 

To show further the role of frequency detuning, we con¬ 
sider the effect on synchronization caused by perturbing 
the adjacency matrix as A^ —Aij (1 + SAuj^). Figure 
[l] reports the synchronization for (a) Erdos-Reyni (ER) 
< k >= 30 and (b) scale-free (SF) networks generated 
by Barabasi-Albert algorithm with < k >= 12. It can be 
seen an enhancement (frustration) of the synchronization 
as S is increased (decreased). Hence, positive (negative) 
values of S potentiate (weaken) the strength of the cou¬ 
plings with AuJij , forcing (preventing) pairs of nodes with 
large detuning to be effectively closer in natural frequen¬ 
cies. Notice that such a local perturbation of the adja¬ 
cency matrix is more effective in promoting or delaying 
the PTS than a global perturbation equally acting on 
all links as shown by the corresponding dashed lines in 
Fig. [lj 

Motivated by Eq. we now merge the structural and 
dynamical information by introducing an effective net¬ 
work characterized by the adjacency matrix C = {C^}, 
where 


@ij — A-ij 


1 - 


A LOi 


Auj n 


(4) 


being Aw maa; the maximum possible detuning present in 
the system. In this way, the effective network exhibits 
the topology of the original network but enhancing those 



FIG. 2. (Color online) Comparison between synchronization 
centrality Ap (black squares) and topological centrality Ap 
(red dots), (a) ER networks, (A;)=50, Ap and Ap are re¬ 
ported vs. the nodes’ natural frequencies c (b) SF networks, 
(A;) = 12, Ap and Ap are plotted vs. the node degrees ki. All 
data refer to ensemble averages over 100 different network 
realizations. 


more synchronizable pairs of nodes according to Eq. 

We remark that, although we are here referring to the 
Kuramoto model of Eq. 0 , Eq. © can be applied much 
broadly to any kind of oscillator’s ensemble which can 
be associated to a set of well defined natural frequencies 
Ui. Indeed, numerical results (not reported here) indi¬ 
cate that the main conclusions we will draw for Eq. 0 
are in fact valid also for networks of chaotic oscillators, 
provided that the power spectrum of each unit is pro¬ 
nouncedly peaked around a unique frequency, i.e. the 
units are strongly phase coherent 22 • 

Now, in order to quantify the role of each node in the 
synchronization process, we extract the most important 
nodes in the effective network defined by C, i.e. we cal¬ 
culate the eigenvector centrality [lj of C, and obtain the 
synchronization centrality vector A c . The i-th compo¬ 
nent Ap > 0 provides a measure of the importance of the 
node i in the effective network and quantifies its potential 
to behave as a seed of synchronization. 

As a testing ground, Fig. [2] shows the comparison be- 
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tween A c and its topological counterpart A a , the eigen¬ 
vector centrality extracted from the original adjacency 
matrix A. For Erdos-Renyi (ER) networks [23] . the 
distribution of the components of the vector A c as a 
function of the corresponding node’s natural frequencies 
shows the existence of many seeds of synchronization 
with natural frequencies close to fl s = 0 (see the black 
squares of Fig. [2ja)). This allows the characterization 
of the connection between the microscale (detuning of 
the links) and the macroscale (emergence of global syn¬ 
chronization) of the system in a much better way than 
A a , whose components (red dots in Fig. [5] a) are instead 
uniformly distributed. For heterogeneous scale-free (SF) 
networks [24] , the synchronization seeds are the hubs and 
therefore A c and A a provide essentially the same infor¬ 
mation, as it can be seen in Fig. [2|b). 

In order to further show how the effective network al¬ 
lows to predict the leading role of each node, we pro¬ 
ceed with confronting our approach with a dynamical 
exploration of the system for the case of ER networks. 
Precisely, we calculate the local synchronization matrix 
S = = \(e lA6ij ) t \ (being A Oij = 6{ — 0j) for a range 

of a values around the synchronization threshold. The 
eigenvector centrality of S, denoted by A s , provides the 
actual dynamical centrality of each node in the synchro¬ 
nization process. In Fig.[3]we report the percentage of co¬ 
incidence between the third of the nodes with the highest 
(lowest) dynamical and synchronization centralities (A f 
and Ap respectively), and the percentage of coincidence 
between the dynamical and topological (Ap) centralities. 
It can be seen that the ranking based on Ap (curves with 
solid symbols) is able to predict up to 65% (75%) of the 
nodes with the highest (lowest) synchronization central¬ 
ity, while the topological centrality (curves with empty 
symbols) only detects at most the 30%. The maximum 
of predictability is reached around the synchronization 
threshold and decreases rapidly due to the homogeniza¬ 
tion of the matrix S for overcritical couplings. 

As synchronization centrality reveals itself as a suit¬ 
able measure to characterize the PTS, we move on elu¬ 
cidating how this quantity helps us also to understand 
the microscopical mechanisms underlying ES. As it has 
been remarked in the introduction, ES can be induced 
when topology and dynamics are related in several spe¬ 
cific ways [TTTIT51 H71T2Q1 [25]. Almost all these meth¬ 
ods are based on a manipulation of the adjacency matrix 
and/or the links weights, such that in Eq. (Jl| A^ is re¬ 
placed by a certain matrix Qjj which usually correlates 
the structural and dynamical features of the network. 

To show how the different studied procedures impact 
the synchronization centrality vector, we compare the 
original A c obtained from Eq. Q with the A^ result¬ 
ing from the corresponding modified effective network 
Cij = fb 7 -(l — ). Results for some of the different 

methods are condensed in Fig. PJ for ER (left panels) and 



FIG. 3. (Color online). Average percentage of coincidence 
between the third of the nodes with the highest (red) and low¬ 
est (blue) dynamical Af and synchronization Ap centralities 
(solid symbols) or topological Ap centralities (empty sym¬ 
bols) . Calculations are performed on ER networks with same 
specification as in the Caption of Fig. [2j and refer to 20 real¬ 
izations of the network’s topology and frequency distribution 
(see text for details). 

SF (right panels) networks. In each panel, Ap (red dots) 
are plotted together with Ap (cyan dots) to show how 
the explosive method actually modified the synchroniza¬ 
tion centrality vector and, therefore, the dynamical role 
of each node. Nodes are sorted in ascending order of Ap. 
In all the cases where the structural and dynamical cor¬ 
relations introduced through Qij successfully lead to ES 
(Figs. [4|a-d,f)), there is an increase (decrease) of Ap of 
those nodes whose Ap was low (high), that is, the weight¬ 
ing method produces a flattening of Ap. In this way, the 
potential ability of the nodes to behave as seeds of lo¬ 
cal synchronization is frustrated until a certain coupling 
strength is reached. Only once the coupling strength is 
large enough, the rest of the network fulfills the condition 
<©, and therefore a sudden transition to synchronization 
takes place. 

More in detail, Figs. |4ja-b) correspond to the method 
described in Ref. [18], where ES is achieved choosing 
Qij = Aij\uJi — ujj | for ER networks (a) and Qij = 
Aij\cji — Uj\hj/ hj f° r SF networks (b), being lij the 
edge betweenness [T]. Figures |4^c-d) show, instead, the 
case Vtij = Aij\(Ji\/ki proposed in Ref. [25] for uniform 
frequency distributions centred in zero. It is easy to see 
that the above increase-decrease compensation is fulfilled 
for both ER (c) and SF (d) networks. As expected, in 
SF networks the modification affects mainly the hubs, 
thwarting their dynamical influence as seeds, and frus¬ 
trating the PTS. Finally, Fig. [4](f ) reports the case of ES 
induced in SF networks by imposing a frequency-degree 
correlation uji = ki [12]. Here the effect is focused on 
the hubs (see the inset), whose synchronization central¬ 
ity is now strongly decreased, while the imposed corre¬ 
lation does not produce a substantial difference between 
A c and A c for the rest of the nodes. There are how¬ 
ever cases when, even if the structure and dynamics are 
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FIG. 4. (Color online). Synchronization centrality of explo¬ 
sive networks (Af, red clouds) versus that of non-explosive 
networks (Ap, cyan clouds) for ER (left panels) and SF (right 
panels) topologies (see text for definitions), (a) and (b) ac¬ 
count for the weighted method of Ref. [18] with a uniform 
natural frequency distribution in [—0.5, 0.5]. (c) and (d) cor¬ 
respond to the weighted method of Ref. [25] with the same 
uniform frequency distribution (see text for the methods’ de¬ 
scription) . Panel (e) is the same as in (d) but natural frequen¬ 
cies are here chosen from the positive values of a Gaussian 
distribution, in a way that ES is no longer induced. Finally, 
(f) corresponds to a SF network where the degree-frequency 
correlation method from Ref. m is applied and the same uni¬ 
form frequency distribution as for panels (a-d) is used. The 
inset in (f) is a zoom centered on the highest degree nodes. 
All data refer to averages over 100 realizations. 

correlated, ES does not occur. For instance, Fig. ge) re- 
ports the same case presented in Fig. [4](c) but for positive 
definite frequencies. And indeed, for this particular fre¬ 
quency distribution, it is seen that the weighting method 
fails to flatten sufficiently the synchronization centrality 
(the red horseshoe cloud), with the consequence that ES 
fails to emerge as well. 

An important application of our measure is the possi¬ 
bility of engineering an efficient strategy to produce ES 
in a generic network by only acting upon a small fraction 


of its nodes, according to a given ranking defining their 
role as synchronization seeds. We here test four possible 
rankings: i) the synchronization ranking based on Ap, ii) 
the distance ranking , that sorts the nodes according to 
the distance AAp = |Ap — Ap|, in) the topology ranking , 
based on Ap and, finally, iv) a random ranking, which is 
used for comparison. These specific rankings are actually 
suggested by the characteristic PTS occurring in both ER 
and SF networks, where the increase-decrease condition 
and the dominant role of hubs constitutes, respectively, 
the essential feature (see Fig. |4j). 

Figure [5] reports the largest jump in the synchroniza¬ 
tion curve r(a) when several percentages of nodes (prop¬ 
erly selected according to the ranking specified in the 
legend) are affected by the weighting methods of Refs, 
nans]. For homogeneous ER networks (Fig. i a)), we 
choose the method from Ref. [THEE]. The synchroniza¬ 
tion ranking (black squares) indicates that a significant 
explosive effect is obtained in the network already for just 
6% of the nodes, whereas the distance ranking (red cir¬ 
cles) requires up to 15% to get an equivalent jump, induc¬ 
ing a complete explosive transition only for percentages 
larger than 30%. In comparison, using a random rank¬ 
ing (blue stars), it is necessary to manipulate at least 
the 40% of the nodes. For heterogeneous SF networks 
(Fig.§b)), we apply instead the degree-frequency corre¬ 
lation method d, choosing a percentage of nodes and 
setting their natural frequencies equal to their degree. In 
this case, both the synchronization (black squares) and 
topological (magenta diamonds) rankings clearly outper¬ 
forms the distance ranking (red circles) by only affecting 
the 10% of the nodes, while the distance ranking behaves 
slightly worst for large percentages. In comparison, the 
random ranking is not able to induce ES even for per¬ 
centages above the 40% d- The differences between 
the two cases are mainly due to the different ways the 
seeds are spread in the network. In the homogeneous 
case, the synchronization ranking obtains an optimum 
effect for small percentages since it focuses on the seeds 
of synchronization. As soon as this percentage increases, 
the nodes with lowest synchronization centrality are not 
captured by such a ranking, and the increase-decrease 
condition is not fulfilled. As there are multiple randomly 
distributed seeds, the distance ranking is only slightly 
better than the random one, as both satisfy the increase- 
decrease condition once the percentage is large enough. 
In the heterogeneous case the seeds are just a few hubs, 
allowing to induce ES acting upon a very small fraction of 
the nodes of the network, whereas the random targeting 
is definitely not the suitable choice. 

In conclusion, we have introduced an effective network 
whose topological properties characterize quantitatively 
the PTS of networked oscillators, connecting the micro¬ 
scopic and the macroscopic behavior of the system dur¬ 
ing the process of synchronization. With our approach 
we can reveal the inner mechanisms beneath ES, which 
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FIG. 5. (Color online). Maximum jump size in the synchro¬ 
nization curve cr(r) vs. fraction of perturbed nodes chosen 
along several rankings: synchronization (back ■), distance 
(red •), topology (magenta ♦) and random (blue *). (a) ER 
networks ( k) = 30, weighting method from Ref. [18] . (b) SF 
networks (k) = 6 imposing uji — ki for the selected nodes 
as in Ref. m- In all cases, data refer to averages over 50 
realizations. 


is shown to be rooted in a frustration of the PTS. This 
approach also provides a very simple yet effective way to 
predict whether or not a dynamical network synchronizes 
explosively for most of the various weighting procedures 
so far considered. Finally, synchronization centrality al¬ 
lows to induce such behavior locally, since we have the 
means to identify and isolate those seeds involved in the 
emergence of synchronization. 
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